Modeling of optical aberration correction using a liquid crystal device 


Xinghua Wang 3 , Bin Wang 3 , ,Paul F. McManamon b ,John J. Pouch 3 , Felix A. Miranda 3 , 

Philip J. Bos 3 * 

3 Liquid Crystal Institute, Kent State University, Kent, OH 44242 
b Air Force Research Lab, Dayton OH 45424-4653 
C NASA Glenn Research Center, Cleveland, OH 44135 


ABSTRACT 

1 3 

Gruneisen ' , has shown that small, light weight, liquid crystal based devices can correct for the 
optical distortion caused by an imperfect primary mirror in a telescope and has discussed the efficiency of 
this correction. In this paper we expand on that work and propose a semi- analytical approach for 
quantifying the efficiency of a liquid crystal based wavefront corrector for this application. 


1. INTRODUCTION 


Researchers at the Air Force Research Laboratory and other groups have pioneered in developing 
liquid crystal devices for tip-tilt, and more general wavefront correction 14 ' 1 1] . McManamon and Watson et 
al. pioneered in liquid crystal Optical Phased Array device, as well as the basic efficiency theory related to 
such device 141 . Gruneisen et al, pioneered in high-resolution wavefront control in a large aperture telescope 
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with diffractive optical element. They studied the diffraction efficiency, dynamic range of correction. 


efficiency for correcting large defocus and astigmatism, wavelength dependence of such system using 
Fourier Optics approach^ l_3 l Also, Gruneisen has addressed the efficiency issue associated with the local 
wavefront slope needed for large correction. 

Here, we present a simple model of evaluating the perfonnance of a liquid crystal based wavefront 
corrector for the telescope correction application. Independent efficiency factors are introduced and 
discussed. 

2. A SIMPLE ANALYTICAL APPROACH FOR ESTIMATING THE CORRECTION EFFICIENCY 


We consider four major factors that contribute to the strehl ratio of the system after correction, as shown 
in Equation 1 . 


1 total ~ 7 7 s 1 ! ray 7 ! ' aber 7 ) ' SLM 


( 1 ) 


The first factor, 17 s , is associated with the passive loss in the system such as the imperfection in the 
mirror surface, scattering, absorption and front surface reflection loss etc. This factor mainly depends on 
the material and process that is used to make the optical elements. 

The second efficiency factor ,rj , is related to the non-conjugate position of the element being 
corrected relative to the corrector. In an ideal perfect conjugate system, where the correction device is 


2 



placed at the conjugate position of the primary optical element, the position of any ray that hits a particular 


point on the primary mirror will hit a corresponding point on the LC-SLM. There is a one to one 
relationship of the ray interception position on the primary mirror and the correction device, and it doesn’t 
matter what is the magnitude or order of the aberration on the primary mirror. In this case, rj my =1 always 

holds true. However, this ideal conjugate relationship may not always hold true for a high numerical 
aperture telescope with large magnification. In this case, the wavefront aberration could cause substantial 
shift of the ray interception position on the correction device for different aberrations, and the efficiency 
rj ray: will be low. While other tenns that govern the overall efficiency increase as the resolution of the 

correction device increases, this tenn drops because the significance of the non-conjugate nature of the 
system is greater if the considered resolution is higher. In can be seen that a ray will be more likely to shift 
to a wrong pixel on the correction device, if the pixel size of the correction device is small; which means, 
this efficiency term drops as the resolution of the correction device increases. 

Not having an approach to quantify the efficiency tenn rj my with in a completely analytical method, 
we can use a simple ray tracing algorithm to describe this term in a physically meaningful way. Consider a 
certain aberration in a certain optical system, where we can trace a mxm bundle of rays throught it. The 
correction device in the system has mxm pixels across the clear aperture. For each pixel on the correction 
device, we can count how may rays intercepts the correction device on that pixel. When there is no 
aberration present in the system, the number of rays intercept a pixel on the correction device should be 
one. However, when aberration is present in the system, the number of rays intercepting a certain pixel is 
no longer one. We define n 0 to be the number of pixels on the correction device that have no ray 
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intercepting that pixel; rp to the the number of pixels on the correction device that have one ray 
intercepting that pixel; n 2 , to be the number of pixels on the correction device that have two rays 
intercepting that pixel; and nk, to be the number of pixels that have k rays intercepting that pixel. Since the 
efficiency of the system will be a direct function of how many rays will miss the desired pixel, which is 
simply ju Q , it can be expressed in Equation (2). 

Iray = 1 “ A() ( 2 ) 


The third efficiency factor//^, is related to the resolution of the corrector. This term becomes larger as 
the resolution of the device becomes higher. If we assume an ideal piston-phase-plate correction device 
with finite resolution, the correction efficiency for random high order aberration is not 100%. Consider the 
simplest case of correcting a simple tilt, and assume q is the number of steps for every wave of tilt, the 
efficiency factor rj aber is the same as the diffraction efficiency for a stair-like blazed grating as in Equation 
3 [4] . 


laberr ~ 0 


si n(;r / q) y 
n / q 


( 3 ) 


We would like to estimate the efficiency factor that, like Equation 3, is related to the resolution of the 
correction device, but is applicable to the more general case of an arbitrary aberration.. Any aberration on 
the primary mirror will induce a corresponding but different aberration on the correction device. The 
aberration profile on the correcting device can be obtained by diffractive beam propagation of a Guassian 

beam from the primary optics to the correction device.. The whole aperture is divided into lxl small 
regions. Within region Cl, if the region is very small, the aberration p^ x y ) ^ can be approximated 
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by a simple tilt, where the maximum difference in the phase across the region is: 

A Q =ma x(P(x,y) Xy y G n )-min(P(x,y) xyeQ ) (4) 

The number of steps along either horizontal or vertical step is S=m/1. Then the number of steps for one 
wave of tilt (q) in region Q is approximately 

Tfl 

<7q=SMq=— (5)' 

/a q 

With such approximation, the efficiency tenn i) aherr , n within region Q can be estimated for a 2-D 
aberration. 


sin (x/q n ) 2 ( sm(xlA a /m) 2 

T laberr,Q. V / z v i \ / ' 

x/q n nlA a / m 


( 6 ) 


If we assume the fractional intensity of the beam on region Q is Wq, then the total efficiency across 
the whole aperture is the weighted average of tj aberr , n across the whole aperture as in Equation (7) [12] . 


7*, - -Z (7) 


xlq <: 


7zlA n / 111 


it is also clear that if the resolution of the correction device is high enough, this efficiency ij aber can 


be very high. 


The fourth factor ij SLM is the efficiency associated with the non-ideal nature of the correction device, 
the LC SLM 123 ' 14 ' 15 l Consider a simple case of a tip-tilt corrector, to correct for hundreds of waves of 
aberration, the thickness of the LC-SLM may become impractically large if the LC-SLM needs to produce 
such large Optical Path Difference (OPD), because the birefringence of the liquid crystal material is 
limited. For monochromatic light, this limitation is resolved by considering an approach where the phase 
profile generated by the LC-SLM is a modulo 2jt version of the desired phase profile. In this case, the 
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maximum OPD change required at any point on the LC-SLM is only A, (X here is the wavelength of light.). 
The phase profile on a LC-SLM is a stair-like blazed grating with phase reset of 2n. As has been discussed 
extensively by [4, 14 \ most pixels in such stair-like blazed grating can produce a phase profile very close to 
the desired phase profile. However, the orientation of the liquid crystal director cannot undergo abrupt 
change to produce the phase discontinuity at the phase reset region. There is a fly-back region where the 
phase slope is in the opposite direction to the desired phase slope. The efficiency of the LC-SLM r/ SLM is 
related to the ratio of the width of the fly-back region and the grating period as expressed in Equation 8 [3 l 
A power of two in Equation 4 is to take into account the efficiency for coherent light. 



( 8 ) 


Here A is the width of each grating segment, A F is the width of the fly-back region. In a 2-D 
LC-SLM discussed in this article, the exact value of this efficiency term depends on many factors such as 
the cell thickness, liquid crystal material birefringence, aperture ratio of the electrode etc. 


However we find A F is primarily related to the strength of the fringing electric field, which is related to 
the cell thickness of the device. Assume the width of the fly-back region is related to the cell thickness by 
an equation of the form shown in Equation 9: 

A r =“(^Y (9) 

Here a is a constant related to the strength of the fringing electric field, d is the thickness of the liquid 
crystal device, PS is the pixel spacing of the LC-SLM. We can fit the diffraction efficiency data obtained 
by director simulation and LDTD simulation 1 16 using Equation 10. 


6 



( 10 ) 


a(d!PS) P Y 


Here N is the average number of pixels for every wave of aberration. The fitted efficiency agrees 
excellently with the FDTD simulation data, for cases where voltage optimization is carried out, with the 
coefficient a as 0.85, and /3 as 0.7 . 

Finnaly we have: 



*1 total ~ 'H sH 'rayH 'aber 7 ) ' SLM 


= ^7, (1 - w 0 ) 


2X( 


n 


sin(;r/A Q / m ) 

7rlA a / m 


2 / a(d/PsY T 

} V 


(ID 


With the above discussion, it is possible to relate the total optical throughput of the system to many design 
parameters as in Equation 1 1 . For example, the total optical throughput is a function of the scattering and 
reflection loss in the system ( // v ), the extent of ray missing corresponding pixel ( w 0 ), the resolution of the 
correction device (m), the complexity of the wavefront (A n ), the pixel spacing of the correction device 
(PS) and cell thickness of the device (d), etc. 
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